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Abstract. We introduce a local algorithm for Khovanov Homology computations — that 
is, we explain how it is possible to "cancel" terms in the Khovanov complex associated with 
a ("local") tangle, hence canceling the many associated "global" terms in one swoosh early 
on. This leads to a dramatic improvement in computational efficiency. Thus our program 
can rapidly compute certain Khovanov homology groups that otherwise would have taken 
centuries to evaluate. 
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1. Introduction 

The "divide and conquer" approach to computation, as ap- 
plied to knot theory, goes roughly as follows. Suppose a certain 
knot invariant takes an exponential amount of time to compute, 
so on a knot with n crossings {n = 14 on the right), it takes 
roughly C" operations, where Ci is some constant. Suppose 
also that that same knot invariant can also be computed "in 
halves"; i.e., it makes sense to compute a "partial" invariant of 
the left half of the knot, consisting of just n/2 crossings, and 
likewise for the right half of the knot. Then, with some luck, it 
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Figure 1. Delooping. 



takes just operations to do each half and if the assembly of the half-computations into 
the full one is cheap then the whole computation takes IC"^^"^ operations, a lot less than the 
original operations. Of course, by iterating this procedure one can save even more and 
carry out the computation in 4C"''^ operations, or even just SCf''* operations. At the limit 
the computation time becomes linear linear in n, at least if one ignores the costs of cutting 
and of assembly. 

In reality, cutting is indeed cheap but assembly isn't. Often the invariant of each half knot 
has to be quite complicated in order to allow for its pairing with every conceivable "other 
half". This often means that each half-knot invariant must take value in a space whose 
dimension grows exponentially in the number of strands h that connect the two halves (6 = 4 
in Figure ??). Thus each half computation takes at least C\ operations and if luck strikes, it 
doesn't take much longer. Typically we can expect h to be roughly the "width" of the knot 
and since we are in the plane, we can expect h to be around \/n. Thus realistically divide 

and conquer may reduce C" to C/". The latter is still very big, but it is a lot smaller than 
the former. 

In fact, the advantage of "divide and conquer" is so big that it is worthwhile to try this 
approach even if the assembly cost is more than C\ or even if no good estimates for the 
assembly cost at all exist, as is the case at hand in this paper. 

This paper applies the "divide and conquer" approach to the computation of Khovanov 
homology [Kh, BNl]. We start in Section 3 with a quick review of the local Khovanov 
theory of [BN2], which amounts to a definition of "Khovanov homology" for half-knots (i.e., 
for tangles), along with the "horizontal composition" technique necessary for the assembly 
of the invariants of two tangles into the invariant of their side-by-side composition. 

In Section 4 we introduce two simple tools, delooping and Gaussian elimination, that allow 
us to "simplify" the invariants of tangles. These tools are the keys to the whole paper, as 
they reduce the complexity of the Khovanov complex associated with a tangle and thus allow 
for much easier horizontal composition assembly. For the impatient reader, delooping and 
Gaussian elimination are depicted in Figures 1 and 2 here. If you understand these sketches 
you've understood the whole paper. 

In Section 5 we describe our algorithm. As an illustration, in Section 6 we "run" the 
algorithm on the figure eight knot. Then in Section 7 we describe in even faster variant of 
the algorithm and in Section 8 we mention the two available implementations and exhibit 
a sample computation. The final Section 9 quickly explains how the tools used in this 
paper lead to a completely automated proof of the invariance of Khovanov homology under 
Reidemeister moves. 
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Figure 2. Gaussian elimination. 



2. History and Acknowledgement 

When I did the research for my paper [BN2] part of my motivation was to set up the 
framework allowing for "divide and conquer" computations as described in Section 1. But 
by the time I finished writing [BN2] I completely forgot about that part of my motivation. 
I wish to thank Marco Mackaay for reminding me! Further thanks to Jeremy Green for 
implementing the algorithm described here and to Louis Leung and Gad Naot for some 
comments and suggestions. 



3. A QUICK REVIEW OF THE LOCAL KHOVANOV THEORY 



(n+,n_) 



Let us briefly recall the definition of the 
Khovanov complex for tangles, following 
[BN2]. Given an n-crossing tangle T with 
boundary 9T, such as the 2-crossing tangle 
displayed on the right, one constructs an 
n-dimensional "cube" of smoothings and 
cobordisms between them (as illustrated 
on the right). This cube is then "flattened" 
to a "formal complex" [T] in the additive 
category Cobl/i{dT) whose objects are for- 
mally graded smoothings with boundary 
dT and whose morphisms are formal lin- 
ear combinations of cobordisms whose tops and bottoms are smoothings and whose side 
boundaries are / x dT, modulo some local relations. 

The Khovanov complex Kh{T) of T is obtained from [T] by some minor further degree and 
height shifts depending only of the numbers n± of over- and under-crossings in T (see [BN2, 
Definition 6.4]). It is a member in the category Kom(Mat(Co6^/^(9T))) of complexes of 
formal direct sums of objects in Coh^n^dT) and it is invariant up to (formal) homotopies. 

For simplicity we are using as the basis to our story one of the simpler cobordism categories 
Cohl/i that appear in [BN2], rather than the most general one, Cob^n (a fuller treatment 
appears in [Na]). It is worthwhile to repeat here the local relations that appear in the 
definition oiCohl/i (see [BN2, Section 11.2]): 
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Also recall from [BN2, Section 5] that [■], and hence Kh{-), are planar algebra morphisms. 
That is, if Ti and T2 are tangles and D{Ti, T2) denotes one of their side-by-side compositions 
(a side by side placement of Ti and T2 while joining some of their ends in a certain way 
prescribed by a planar arc diagram D), then [D{Ti,T2)] = i5([Ti] , [T2]). Here as in [BN2, 
Section 5] D([Ti],[T2]) is the "tensor product" operation induced on formal complexes by 
the horizontal composition operation D on the canopoly Coh^ji. In exactly the same sense 
we also have that Kh{D{Ti,T2)) = D{Kh{Ti), Kh{T2)). 

Thus Khovanov homology is ready for a divide and conquer computation. It makes sense 
for half-knots (tangles) and there is a composition rule that takes the invariants of the halves 
and produces the invariant of the whole. But as it stands there is no advantage (yet) for this 
computation method. If Tj (for i = 1,2) has crossings, the Khovanov cube for Tj consists 
of 2"* vertices, thus [Tj involves 2"' objects and thus d([Ti] , [T2]) involves 2"^ ■ 2"^ = 2"i+"2 
objects, exactly as many as in [D(Ti, T2)], and nothing has been saved. 



4. The tools: delooping and Gaussian elimination 

To overcome the difficulty from the previous paragraph we need to learn how to simplify 
[Ti] and [T2] (modulo homotopy) before taking their tensor product and thus before the 
biggest number (2"i+"'2^ jg encountered. For this we need tools for simplifying complexes 
over the category Cobl^i. These tools are Lemma 4.1 and Lemma 4.2 below. 

Lemma 4.1. (Delooping) If an object S in Cobl^i contains a closed loop i, then it is isomor- 
phic (in Ma,t{Cobl^i)) to the direct sum of two copies 5"{+l} and S"{ — 1} of S in which i is 
removed, one taken with a degree shift of +1 and one with a degree shift of —1. Symbolically, 
this reads Q = 0{+l} © 0{-l}. 

Proof. Here are the isomorphisms: 




It is easy to verify using (all!) the relations in (1) that the two possible compositions of the 
morphisms above are both equal to the identity morphisms of the relevant objects. (This is 
the only place in this paper where the relations in (1) are used). □ 
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Lemma 4.2. (Gaussian elimination, made abstract) //</>: 61 ^ 62 is an isomorphism (in 
some additive category C), then the four term complex segment in Mat(C) 
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is isomorphic to the (direct sum) complex segment 
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Both these complexes are homotopy equivalent to the (simpler) complex segment 



(5) 
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Here C, D, E and F are arbitrary columns of objects in C and all Greek letters (other 
than (f) ) represent arbitrary matrices of morphisms in C (having the appropriate dimensions, 
domains and ranges); all matrices appearing in these complexes are block-matrices with blocks 
as specified, bi and 62 o'^e billed here as individual objects of C, but they can equally well be 
taken to be columns of objects provided (the morphism matrix) cj) remains invertible. 

Proof. The two 2x2 (block) matrices in (3) and (4) differ by invertible row and 
column operations (i.e., by a "change of basis"). When the corresponding column and row 

operations are performed on (/i v) and on , the results are (/i — h"y(()~^ = (O 

respectively (note that /i0 — z/7 = and ipa — 5/5 = as in 



and 



a 



the original complex the differential squares to 0). Hence the complexes in (3) and (4) are 
isomorphic. 

The complexes in (4) and (5) differ by the removal of a contractible direct summand 

^ bi — ^ 62 ^ (remember that is invertible). Hence they are homotopy equiv- 
alent. □ 



5. The algorithm 

We thus have a procedure for simplifying complexes Q with objects in Cob,/i (and thus, 
along with "divide and conquer", we have a potentially fast way of computing Khovanov 
homologies) : 

• Whenever an object in Q contains a closed loop, double it up using Lemma 4.1, 
removing the closed loop and inserting ±1 degree shifts. This done, Q becomes Q' . 

• Q' is "bigger" than Q, but it is made up of many fewer possible objects. Thus it is 
likely that many morphisms in Q' are isomorphisms. Whenever you find one, cancel 
it using Lemma 4.2. Call the result of doing this iteratively Q". 

A priori, we cannot guarantee that Q" will be simpler than Q. But experimentation shows 
that it is, as seen in the next few sections. 
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6. The figure eight knot 



Our first example is the figure eight knot, cut in half into two tangles 
Tl and Tr as shown on the right. Tl is the tangle XX- Its Khovanov 
complex is the complex Qi = Kh{yC\) appearing below. Here and later 
we use the same notational conventions as in [BN2]. Thus the height zero 
object in the complex is underlined, and X is the saddle cobordism whose 
domain is X and whose range is ) (: 




{-4}- 



{-3} 



X) 



(X 



{-2} 



Only one of the objects in Qi contains a loop — the last one. Delooping it using Lemma 4.1 
we get the complex Q2, which is isomorphic to Qi (we've also replaced some of the smoothings 
and cobordisms appearing in fli with isotopic ones): 



^9 : 



{-4}- 



d- 



)( 



{-3}- 



>( -X 



{-3} 
{-1} 



The symbol ) ( here, when appearing as a cobordism, denotes the identity automorphism 
of the smoothing ) ( . Likewise ) ( and ) ( denote that same cobordism, with an extra dot on 
the left (or right) "curtain". 

The upper left entry in the differential d~^ in Q2 is an isomorphism (so is the upper right 
entry, but one is enough). So we can apply the first part of Lemma 4.2 taking to be 
that upper left entry. The result is the complex below, which is isomorphic to ^2 and 
fli. Note that in this case we have no "(/i i^)" term, and that in this case, — 70~^(5 = 
->(o)(o(-)() =+)(: 







{-4}- 



{-3}- 







-)< + >( 



{-3} 
{-1} 



Dropping the contractible summand as in the second part of Lemma 4.2 we get the complex 
^4: 
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The complex ^l^ contains fewer objects than fii, hence from a computational perspective 
it is indeed simpler. The savings, about 25%, may not appear to be much, but remember — 

• This 25% will be compounded with whatever savings we will incur later in the com- 
putation. 

• This is just a sample. In real life the savings are greater. 

In a similar manner we can compute and simplify the Khovanov complex of the right half 
of the Figure Eight knot, the tangle Tr. The result is the complex ^^4: 




)( {4} 



Next we have to take the "tensor product" of ^4 with \&4 using the same side-by-side 
composition used to make the Figure Eight knot out of Ti and Tr. The result is the double 
complex Ti below. To prevent clutter we've done away with the vector and matrix brackets; 
also, notice the signs in the middle row, sprinkled as commonly practiced on tensor products. 




Ti may seem monstrous, but yet, it contains just 9 objects, as compared to the 16 objects 
one sees in a direct computation of the Khovanov complex of a 4-crossing knot. Anyway, 
it is best to re-write Ti a bit before proceeding. The smoothings and cobordisms can be 
smoothed out, and dots can be moved around cobordisms so as to cancel the four differences 
on the lower left of the Ti diagram. The result is T2 below. 
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If we were a computer program we would have now "flattened" the double complex T2 

to a single complex of the schematic form (■) — >■ (:) — > (:) (:) (■). Such "single" 
complexes are more easily manipulated on a computer. But it is unlikely this paper will ever 
be appreciated by anything but humans. So we'll stick to the more readable double complex 
form while remembering that we really have just one differential going south and east. 

Anyway, the next step is to replace every loop in every object in T2 with a pair of (degree- 
shifted) empty sets, as in Lemma 4.1, while replacing the differentials with their compositions 
with the explicit isomorphisms of (2). But there's nothing but loops in T2, so we are left with 
a complex T3 in which all the objects are degree-shifted empty sets and all the morphisms 
are (matrices of) scalar multiples of the empty cobordism (note that modulo the relations 
in (1) all closed surfaces reduce to scalars): 
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Note in passing that the matrices m and A appearing here are the matrices representing 
the product and the co-product of [Kh] relative to the basis {X, 1) used there (or the basis 
{v-,v+) used in [BNl]). This is essentially the content of [BN2, Section 9.1]. 



We can now apply Lemma 4.2 repeatedly to T3, 
until no invertible entries remain in any of the matri- 
ces. Over Q any non-zero number is invertible so our 
process stops when all matrices are 0. Thus a lengthy 
iterated application of Lemma 4.2 stops at the com- 
plex T4 shown on the right (a human could save quite 
a lot by being clever, but that's not our point here). 

Belatedly flattening the double complex T4 we ar- 
rive at our final answer, the complex 



[0{-5}] — [0{-i}] — 



■0{-iy 

0{i} 



[0{i}] 



[0{5}] 



(6) 



[0{-5}] — [0{-i}] 



■0{-iy 

0{i} 



[0{i}] — [0{5}] 



To recover homology groups out of T5 we need to apply to it some functor JF taking Coh^n 
to graded vector spaces, and then take homology. The latter step (computing the homology) 
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is the do nothing operation as all differentials are 0. Typically (e.g., as in [BN2, Section 
11.2]) the functor T maps the empty smoothing to the one dimensional vector space (call 
it Q) at degree 0. And so we can read directly from (6) that the "conventional" Khovanov 
homology over Q of the figure eight knot is 6-dimensional with generators at bidegrees 
(-2,-5), (-1,-1), (0,-1), (0,1), (1,1) and (2,5). 



7. A FASTER ALGORITHM 



It turns out that there is a somewhat better way to 
assemble our two main tools into a running algorithm. 
Instead of computing each "half knot" and combining 
the result, "scan" the knot from left to right (or top to 
bottom, as in the picture on the left), adding one crossing 
at a time. After each crossing is added, use delooping 
(Lemma 4.1) and Gaussian elimination (Lemma 4.2) to 
simplify the result. 

Since we have no rigorous estimate, we have to settle 
with a possibly naive estimation of the complexity of the 
"divide an conquer" algorithm and of the new "scanning" algorithm. In both cases the 
bottleneck ought to be where the knot is widest; if the width at the widest cut is W , we 
expect the complex corresponding to either half of the knot to be of size for some C . In 
the "divide an conquer" algorithm we need to "tensor multiply" two such complexes so for a 
little while we have to hold a double complex of total size of about (C*^)^ = (7^^ (compare 
with the complex Ti above). In the "scanning" algorithm we only need to tensor multiply 
the complex for the left half-knot with the complex of a single crossing, whose size is 2. So 
we only see a double complex of size 2C^ before we get the chance to simplify again. 

8. Computer programs 




Two implementations of the "scanning" algorithm are available, 
both as a part of the knot theory package KnotTheory' [KT]. The 
first one, FastKh, was written by the author just to test the prin- 
ciple, with no attempt at optimization. Yet for example, it was 
able to compute the Khovanov homology of the 35-crossing (7, 6) 
torus knot T7.6 in about one day of work; a task that would have 
taken about a 1,000 years without the use of tangles. The second 
one, JavaKh was written by Jeremy Green, a summer student of 
the author's at the University of Toronto, in the summer of 2005. 
The computation leading to Table 1, for example, of the Khovanov homology of the 48- 
crossing (8,7) torus knot Tgj, takes a few minutes using JavaKh. The (r, j) entry of that 
table contains the degree 2r -|- j piece of the rth integral Khovanov homology Q2r+jH^(Tsj) 
of Ts 7. 




9. The Reidemeister moves 



One further advantage of the ability to compute with tangles is that the proof of up-to- 
homotopy invariance of the Khovanov complex under Reidemeister move becomes mechan- 
ical. All that one has to do is to compute and simplify the complexes corresponding to 
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Table 1. The integral Khovanov homology groups Q2r+jH'''{T8j). 



each side of each Reidemeister move. One gets the same result for both sides of any given 
Reidemeister move and hence invariance is proven^. 

As an example, Figure 3 proves invariance under the second Reidemeister move i?2 in 
summary form. Start on the upper left, where the "hard side" of i?2 is presented as the side- 
by-side product of two single-crossing tangles. Move on to the right where the corresponding 
double complex is shown, and back to the bottom left where that complex is flattened, 
delooped and the two isomorphisms that appear in the result are canceled out. What remains 
is a single-entry complex equal to the "easy side" of i?2. For simplicity degree shifts are 
ignored in Figure 3. 



In fact, the computer programs discussed above can carry out these computations, turning them hteraUy 
mechanicaL Though the presently available front ends for those programs are only configured to take knots 
as inputs. 
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Figure 3. Proof of invariance under R2. 
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